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Abstract 

An inhomogeneous first-order integer-valued autoregressive (INAR(l)) process is 
investigated, where the autoregressive type coefficient slowly converges to one. It is 
shown that the process converges weakly to a Poisson or a compound Poisson distri- 
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1 Introduction 

A zero start inhomogeneous first order integer-valued autoregressive (INAR(l)) time series 
(X n ) n&+ is defined as 

X n -\ 

X n = E fnj +e„,, neN, 

j=i (1-1) 

X = 0, 

where {£ n ,j, £n '■ n, j G N} are independent non- negative integer- valued random variables 
such that {£ nj - : j G N} are identically distributed and P(£ n ,i £ {0, 1}) = 1 for each 
n G N. In fact, (AT n ) ngZ+ is a special Galton-Watson branching process with immigration 
such that the offspring distributions are Bernoulli distributions. We can interpret X n as 
the size of the n th generation of a population, £ n t j is the number of offspring produced 
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by the j th individual belonging to the (n — l) th generation, and e n is the number of 
immigrants in the n th generation. 

The process (11.11) is called INAR(l) since it may also be written in the form 

X n = g n o X n _i + e n , ne N, 
X = 0, 

where 

Qn '■= E^ n> i 

denotes the mean of the Bernoulli offspring distribution in the n th generation, and we use 
the Steutel and van Ham operator g o which is defined for g G [0, 1] and for a non-negative 
integer-valued random variable X by 

go X : = < j=i 

(o, x = o, 

where the counting sequence consists of independent and identically distributed 

Bernoulli random variables with mean g, independent of X (see Steutel and van Harn 
[19]), and the counting sequences involved in g n o X n _i, nGN, are mutually independent 
and independent of (e n ) n eN- 

Let us denote the factorial moments of the immigration distributions by 

:= Ee n (e n - 1) • • • (e n - k + 1), n, k e N. 
If m n i < oo for all n6N then we have the recursion 

EX n = ^ n EX n _! + m nA , neN, 



since 



'X 



n-l 



3=1 



X 



n-l 



n,j + Ee n — X n _\g n + m n> x. 



3=1 



Consequently, the sequence (g n )n£N of the offspring means plays a crucial role in the 
asymptotic behavior of the sequence (X n ) neZ+ as n — » oo. The INAR(l) process (X n ) ne z + 
is called nearly critical if £> n — > 1 as n — > oo. We will investigate the asymptotic behavior 
of nearly critical INAR(l) processes. 

Non-negative integer-valued time series, known as counting processes, arise in several 
fields of medicine (see, e.g., Cardinal et al. [8] and Franke and Seligmann |13j). To model 
counting processes Al-Osh and Alzaid [3] proposed the INAR(l) model. Ispany et al. 
investigated the asymptotic inference for nearly unstable INAR(l) models. Later on Al-Osh 
and Alzaid [5] and Du and Li [TO] generalized this model by introducing the INAR(p) model. 

The INAR models are special branching processes where the offspring distributions are 
Bernoulli distributions. The theory of branching processes has been developed for a long 



time, see Athreya and Ney [6], and it can be applied in various fields. Branching processes are 
well-known models of binary search trees, see Devroye [9]. A recent application of them is the 
domain of peer-to-peer file sharing networks. Traffic measurements show that the workload 
generated by P2P applications is the dominant part of most of the Internet segments. The 
file population dynamics can be described by these mathematical models which also make 
possible the design and control of peer-to-peer systems, see Adar and Huberman [2], Zhao 
et al. [20]. Space-time processes are standard models in seismology, see Lise and Stella [18] . 
One of these is the Epidemic-type Aftershock Sequence (ETAS ) model and they serve for 
surveillance of infections diseases as well, see Farrington et al. [12]. The theory of branching 
processes can also be applied to data on different aspects of biodiversity or macroevolution 
by the help of using phylogenetic trees, see, e.g., Aldous and Popovic [3] and Haccou and 
Iwasa [16]. An inhomogeneous branching mechanism has been considered in Ispany et al. 
[13] . Drost et al. [H] proved that the limit experiment of a homogeneneous INAR(l) model 
has a Poisson distribution. 

The present paper seems to be the first attempt to deal with the so-called nearly unstable 
inhomogeneous INAR(l) model. The paper is organized as follows. In Section 2 two basic 
lemmas are proved for inhomogeneous INAR(l) process. In Section 3 the case of Bernoulli 
immigrations, in Section 4 the case of non-Bernoulli immigrations with Poisson limit distri- 
bution are considered. Section 5 is devoted to the general case when the limit distribution is 
a compound Poisson distribution. The results are extended for triangular system of mixtures 
of binomial distributions. In the Appendix at the end of paper some technical lemmas are 
gathered. 

2 Preliminaries 

Let Be(p) denote a Bernoulli distribution with mean p £ [0,1]. The distribution of a 
random variable £ will be denoted by £(£)• Consider the unit disk D := {z £ C : \z\ ^ 1} 
of the complex plane C. The (probability) generating function of a non-negative integer- 
valued random variable £ is given by z h- > E(z^) for z £ D, and we have E(zr) £ D for 
all z £ D. Introduce the generating functions 



Lemma 1 For an arbitrary inhomogeneous INAR(l) process {X n ) n& i + we have 



F n (z) := E(z x "), G n (z) := E(z^), 



H n {z) : 



z £ D. 



Fn{z) = Y[H k (l + Q [kM (z - 1)) 



n £ N 



k=l 



for all z £ D, where 




for 1 ^ k ^ n — 1 



for k = n. 



Proof. The basic recursion for the generating functions F n , neN, is 



FJz) 



j=l ?"J +t » 



7 ~r£n 



X 



n-1 



XV, 



=E(G n (z) x ^)H n (z) = F n ^(G n (z))H n (z), 

valid for all z G C with z G D and G n (z) G D, see Athreya and Ney j6j p. 263]. Clearly 
z G D implies G n (z) G D, hence ( 12.1 p is valid for all z E D. Since £(£ n ,i) = Be(f? n ), we 
have 

G n {z) = 1 - £ n + Q n z = 1 + n (z - 1) 

for all z G C. We prove the statement of the lemma by induction. For n = 1, we have 
F\{z) = Hx(z) = Hi(l + (z — 1)). By the recursion (12. ip , we obtain for n ^ 2 

n-1 

F n (z) = F„_i(l + g n (z - l))H n (z) = H n (z) ]jH k (l + g n Q[ k , n -i](z - 1)) 



k=l 



Y]_H k (l + Q[k, r 



D): 



fc=l 



□ 



and the proof is complete. 

In fact, X n can be considered as a sum of independent Galton- Watson processes without 
immigration. Namely, 



IX 



n G N, 



(2.2) 



fc=i 



where 



Yr 







for = 0, 



71, fc 



Y"n-l,l + -+>n-l,* 

£ n j for 1 ^ /c ^ rt — 1, 

j=K„_i,i+-H-y„_i,fc_i+i 

e„ for k = n. 



(2.3) 



The distribution of Y njk is a mixture of binomial distributions with a common probability 
parameter g n , since the number Y n _i ik of Bernoulli random variables in the sum (12.31) is a 
random variable as well. For a probability measure fi on Z + and for a number pG [0,1], 
the mixture Bi(/x,p) of binomial distributions with parameters and p is a probability 
measure on Z + defined by 

Bi(/*,p){j} := (1 -pf- J m for J G Z+. 

1=3 

It is a particular example for mixture of distributions, see Johnson and Kotz [TTJ Section 
1.7.3], because the common method for mixture of binomial distributions is to use differ- 
ent values of probability parameter, see Johnson and Kotz [T71 Section III. 11]. Note that 
Bi(/i, 1) = \i. 

Lemma 2 For all n G N, 1 ^ k ^ n, the distribution of Y njk is a mixture of binomial 
distributions with parameters e k and Q[ k , n ] ■ Thus 

C{X n ) = * Bi Q[ k , n ]), 

where * denotes convolution of probability measures. 



Proof. First we check that Bi(Bi(/x,p), g) = Bi(ju, pq) for an arbitrary probability measure 
/i on Z + and for all p, q G [0, 1]. Indeed, for all j G Z + , 



Bi(Bi(/i,p),g){j} = - g )^Bi(/,,p){l} 



oo 

/ / \ / / - \ 

J/1 „\fc-i 



/(i-pr^} 



oo fc 

\k-i 



f>)>mE f £ ) ( 1 f)( P (i-q)Y- j (i- P y 

k=j e=j v/ \ / 

k=j t=j V ' 



k=j 

= Bi(/i,pg){j}. 

Since Y^ x = e l5 thus £(Yi,i) = Bi(£(ei), 1), and £(F n ,fc) = Bi(£(Y„ ife _i), g n ) for all 
n ^ 2 and all k — 1, . . . , n, we obtain the statement of the lemma by induction using the 
previous argument. □ 

Remark that Lemma [2] implies the formula given for the generating function of X n in 
Lemma [U since the generating function of a distribution Bi(/x,p) is z i— > if (1 + — 
where if denotes the generating function of /i. 

3 Poisson limit distribution: the case of Bernoulli im- 
migrations 

First consider the simplest case, when C(e n ) = Be(m n> x), nGN. 

Theorem 1 Let (X n ) n& % + be an INAR(l) process such that P(e n G {0,1}) = 1 for all 
n G N. Assume that 

oo 

(i) g n < 1 for all neN, lim g n = l, (1 _ 0n) = oo, 



n=l 



ii) lim^ = A6[0,(X)) 



... W 

Then 



X n Po(A) as n — > oo. (3.1) 



(Here and in the sequel Po(0) zs understood as a Dirac measure concentrated at the point 
1-) 



Remark that the condition (1 — g n ) = oo may be replaced by f\ Qn — 0. Moreover, 

n=l n=l 

if A > then the condition lim g n = 1 may be replaced by lim m n i = 0. 

n— *oo n—*oo ' 

Proof. In order to prove the statement, we will show that 

lim F n (z) = e x(z ~ 1] 

for all z G D. Since C(sj) = Be(mj,i), we have 

H k (z) = l + m kt i(z-i), zeD, keN. 
Applying Lemma [TJ we can write 

n 

F n (z) = l[[l+m kil Q [kjri] (z-l)], zeD, neN. (3.2) 
fe=i 

Consider the functions F n : C — ► C, n G N, defined by 

n 

F n (z) -=\\e m ^ e ^ iz - l) . (3.3) 

k=l 

In fact, f!3.3p is the generating function of a Poisson distribution. The terms in the products 
in ( 13.21) and ( 13.31) are generating functions of probability distributions, hence Lemma [3] is 
applicable, and we obtain 

n 

\F n (z) - F n (z)\ *C \e mh ' ie ^ {z - 1] - 1 - m kA Q [ktn] (z - 1)| 

k=l 

for z £ D, nGN. An application of the inequality \e u — 1 — n|^|w| 2 valid for all u G C 
with \u\ ^ 1/2 implies 

WD - 1 - m k>lQ[kM (z - 1) < m^4 )n] |z - 1| 2 (3.4) 



for zGC with wi^ i |# — 1| ^ 1/2. By Lemma [5] and taking into account assumption 

tn n ,i 



lim = A G [0, oo), we have 



max mfc^^jfc^] = max a nk — >• U as n->oo. (3.5) 

Thus, the estimate (13. 4p is valid for all z G D, for sufficiently large n and for all 
k = 1, . . . ,n, and we obtain 



k=l 



\F n (z)-F n (z)\^\z-l\ 2 J2™l,ieh 

Til? 1 

By lim = and by Lemma [5] we obtain 



2 

m 



X] m M^M = ~r~T a ni as n ^ oo. (3.6) 



Consequently, 

lim \F n (z) - F n (z)\ = for all z E D. 

n— +00 

An application of Lemma yields 

n n 

5_TOk,i0[M = ^2 1 _ fc;1 a S A as n ^ 00. (3.7) 

Consequently, 



k=x k=i ~ 0k 



lim F n (z) = e Mz - 1] for all z E D, 

n^oo 

and we obtain F n (z) — > e^ 2-1 -* as n — ► 00 for all 2 G D. □ 

Second proof of Theorem [I] by Poisson approximation. We may prove the theorem 
by Poisson approximation as well. The total variation distance between two probability 
measures \i and v on Z + equals 



2 

i=o 

A sequence (yU ra )neN of probability measures on Z + converges weakly to a probability 
measure /i on Z + if and only if <i(/i n ,/i) — > 0. We prove (13. ip by showing that 

d(£{X n ),Po(X)) -»• as n^oo. (3.8) 

One can easily check that Bi(Be(p), gj = Be(pq) for arbitrary p,q E [0, 1], hence by Lemma 
[2] we obtain £(X n ) = * Be (pik,iQ[k,n]) ■ By Lemma IU 



d\C(X n ), * Po(m fc ,ie[fc )tl ])J < 5^^(Be(m fcj i^[ fej „]),Po(m fe) i^[ fe)ri ])). 

fc=i 

We show that 

d(Be(p),Po(p)) ^p 2 (3.9) 

for all p E [0,1]. Indeed, 

d(Be(p),Po(p)) = ^(e- p -l+p) + ^(p-pe^) + ^(l- e - p -pe^) = p{l - e~ p ) ^ p 2 . 
Applying (I3.9P and (13.61) . we conclude 

n 

d(c(X n ), * = ¥o{m Kl Q [Kn] )} < J]my M] -»• 0. 

Clearly, 



fe=i 



* 

fc= 



Po(m fc) i^[ fc) „]) = Pof y^m fc;1 g[ fc)n ] j -> Po(A) 



______ 23 

in law by (13.71) , and we obtain X n — ► Po(A). □ 



4 Poisson limit distribution: the case of non-Bernoulli 
immigrations 



Theorem 2 Let {X n ) n& i + be an inhomogeneous INAR(l) process. Assume that 



(i) g n < 1 for all neN, lim g n = 1, ^ (1 — g n ) = oo, 



n=l 



ii) ii m = \ e [0, oo), lim = 0. 



Then 



Remark 1 Since 



X n -—»■ Po(A) as n — > oo. 



n,l = J^P(£n>i), m n , 2 = 2^jP(e n 



assumption (ii) implies 



lim P ' £ " > 2 » = 0. 

n^oo 1 — Q n 



(4.1) 



In general the converse is not true. However, if there exists a sequence (&j)jeN of non- 
negative real numbers such that Y^Li jbj < oo and p< f n> ^ ^ ft • for all ?', n G N, then 
(14. ip implies (ii) by the dominated convergence theorem. 

Proof. By Lemma [U we can write 

n 

F n (z) = Y[H k (l + g [k>n] (z-l)), zeD, neN, 

k=l 

Consider the functions F n : C — > C, neN, defined by 

n 

F n {z) = Y[ [1 + m Kl g [Kn] (z - 1)] . 
fc=l 

By Lemma [3j we obtain 

n 

\F n (z) - F n (z) \ < ^ Lff fe (l + ^[ fc , n] (2; - 1)) - 1 - m k ,i Q[ k ,n]{z - I] 



k=l 



for z G -D, neN. Applying Lemma El we have 



1 



Thus 



\H k {u) - 1 - mfe,i(« - 1)| < -m fc , 2 |w - 1| 2 n£D, fceN. 



#fc(l + ?[&,n](^ - 1)) - 1 - £[fc,n]0 - 1) < 2 mfc ' 2 ~ 1 



for all z G D, since z G D implies 1 + Q[k,n]{z — 1) G -D. Consequently, 

1 n 

|F n (z) - F n (*)| < -\z - l\ 2 J2 m ^ 4,n] - 



fc=l 



as n —>■ oo for all z G D by Lemma [5] taking into account assumption lim = 0. 

n— >oo 1 

Theorem [1] clearly implies F n (z) — > e A ^ _1 - ) for all z E D, hence we conclude F n (z) — > e A ^ -1 ) 



as n — > oo for all z & D. 



□ 



Second proof of Theorem [2] by Poisson approximation. Note that m k ,iQ[k,n] ^ 1 

for sufficiently large n and for all 1 ^ k ^ n by (13. 5ft . By Lemmas [2] and HI we have, for 
sufficiently large n G N, 



d\C(X n ), J^Be(m fcjl £[ fcjn] )J < ^ d(Bi(e k , Q[k,n]), Be(m M £[ fc)n ])). 

fc=i 

We prove that 

d(Bi(e,p),Be(pEe)) < ^ 2 Ee( £ - 1), (4.2) 

where p G [0, 1] and £ is a non- negative integer- valued random variable such that pEe ^ 1. 
We have 



d(Bi(e,p),Be(pEe)) ^±(A + fl + C), 



where 



.4 



B 



J2(l- P yP(e = £)-(l-pEe) 

1=0 

oo 

^£p(l - pY~ l P(e = t) - pEe 



oo oo 



< J2\{l-pf-l+ip\ P(e 

oo 

^p^|(l-p) € - l -l|P(e = 



/A /A 

j=2 fc=j £=2 j=2 

oo 



*£p(s = e)(i-(i-py-£ P (i- P y- 1 ) 



1=2 
oo 



^ J] P(e = - £p - (1 - pY\ + fy\l - (1 - p)^). 



£=2 



By Taylor's formula for the function p \— > (1 — we get 

|1 (1 < -£(£- l)p 2 sup (i-^/- 2 <^-iy. 

2 0e[o,i] 2 



|(l-p)^_l|^p(i-l), 



Finally, since 

we obtain (14.21) . Thus, we have 



k=l 



where the right hand side tends to by the assumption lim ^hhl. = o. Obviously, Theorem 
[T] implies 

d ( ^Be(m k)1 g[ k)n ]), Po(A)j ^ as ii^oo, 

hence we obtain 

d(£(X n ),Po(\)) ^ as n^oo, 

which completes the proof. □ 

In fact, a similar theorem holds for triangular system of mixtures of binomial distribu- 
tions. 

Theorem 3 Let k n G N for all jigN, and {C n ,k '■ 1 ^ k ^ k n , n G N} be non-negative 
integer-valued random variables. Moreover, let p njk G [0,1], 1 ^ k ^ k n , n G N. Assume 
that 

(i) E Pn,kE(n,k -> A for some A ^ 0; 
fc=i 

(ii) E (Pn,fcEC„, fc ) 2 -> 0; 
fe=i 

(iii) E Pn. fc EC„, fc (C„, fc - 1) - 
as n — > oo. T/ien 

* Bi(£(C n) fc),Pn,fe) -»• Po(A) 

in law as n — > oo. 



5 Compound Poisson limit distribution 

Recall that if /i is a finite measure on Z + then the compound Poisson distribution CP(fi) 
with intensity measure \x is the probability measure on Z + with generating function 

z i — ► exp < il{j\{z> — 1) > for 2 G D. 

In fact, CP(/x) is an infinitely divisible distribution on Z + with Levy measure /i restricted 
onto N, and, for an arbitrary infinitely divisible distribution v on Z + , there exists a 
finite measure /i on N such that z/ = CP(/x). Moreover, CP(/x) is the distribution of 
the random sum 

n 

£<., 

i=i 

where {II, ^ : j G N} are independent random variables, C(H) = Po(||/x||) and £(£j) — tutt 
for j G N, where := Ejli/^O}- Further, CP(/x) is the distribution of the weakly 
convergent infinite sum 

oo 

52 w> 



where {rjj : j e N} are independent random variables with C(rjj) = Po(fi{j}) for j e N. 
(See Barbour et al. Section 10.4].) 

First we consider the case when the intensity measure fi of the limiting compound 
Poisson distribution CP(/i) has bounded support. 

Theorem 4 Let (X n ) ne z + be an inhomogeneous INAR(l) process. Assume that 

oo 

(i) g n < 1 for all nGN, lim g n = l, ]T (1 - g n ) = oo, 

(ii) lim j^-s = X i e [°> °°) / or j = 1, • • • , ^ Aj = 0. 

n— »oo J 1 . 1 

Then 

X n — ► CP(/i) as n — > oo, 
where \x is a finite measure on {1, . . . , J — 1} given fry 

1 J ~ i_1 f-lV 

: = - E ^Vi. j = 1, . . . , J - 1. (5.1) 

Remark 2 One can easily check that Aj, j = 1, . . . , J — 1, are the first J — 1 factorial 
moments of the measure /i, i.e., 

j-i 

x j = Y^ i ^- l )---^-3 + ^Oh j = i, . . . , j — i. 

Moreover, since 

°° / — 1\ 

assumption (ii) implies 
(ii)' lim = G [0, oo) for z = 1, . . . , J with fx{J} = 0. 

?1— >0O l ^ @ n > 

On the other hand, (ii)' and additional domination assumption, see Remark [TJ imply (ii). 
Proof. By Lemma [TJ we can write 

n 

F n (z)=]jH k (l + Q [ktn] (z-l)), zeD, neN. 

k=l 

Consider the functions 

n 

F n (z) := JJ e Hk(i+e lkM ( z - l »-\ zE D, neN. 
fc=i 

By Lemma [31 we obtain 

n 

\F n (z) - F n (z)\ < ^ -H k (l + g [kM (z - 1)) 

fc=i 



for z G D, nGN. An application of the inequality \e u — 1 — u|^|m| 2 valid for all u G C 
with \u\ ^ 1/2 implies 

e H k (i +C[k , n ,(*-i))-i _ Hk ^ + Q{k ^ z _ 1} ) <^ Hk ( 1 + g M ( z - 1)) - 1 2 (5.2) 



for 2 G C with 



Thus 



-£/fc(l + £>[fc,n](2 — 1)) — 1 ^ 1/2. Applying Lemma El we have 
\H k (u)-l\^m kt i\u — l\, ueD, keN. 

H k (l + Q[ k ,n]{z ~ 1)) - 1 < "lfc,l 0[*,n]l Z ~~ X l 



for all z £ D, since z E D implies 1 + Q\h, n ]( z — 1) G -D. By Lemma [5] and taking into 
account assumption lim = Ai G [0, oo), we obtain (13. 51) . Thus, the estimate (15.21) is 
valid for all z £ D, for sufficiently large n and for all k — 1, . . . , n, and we obtain 

n 

|F n (z) - F n (^)| ^ 1 2; — 1| 2 E m /M ^ 2 fc n] — > as n — > 00 for all z G D 

fc=i 

by dSSD- Clearly 

F„(z) = exp <^ ^ Jif fe (l + e[ fc , n ](2 - 1)) - 1 
I fe=i 

Again by Lemma [6], we have 



j-i 



H k {l + g [k , n] (z -1))-1 = J2 -r4«]( z _ 1)J + ^' fc ' j( 

3=1 J ' 

for all z G D, for sufficiently large n and for all k = 1, . . . ,n, where 

\Rn*j{z)\<^jfl!lk,n)\z-l\ J - 

An application of Lemma [5J yields 

n n 

Em J m fc.3 Jj) _^ X 

fc=l fe=l 

for j — 1, . . . , J. Moreover, by (15.31) . 

n 



as n — > 00 



(5.3) 



fe=l * fc=l 



as 7r — > 00 for all z E D since Aj = 0. Consequently, 



j-i 



.3=1 



lim F n (z) = exp <^ V (z - 1) J J> =: F(z) for all 2 G D, 



where F is the generating function of a probability distribution. Clearly 

J-l > J-l , 3 / -\ J-i i J -i 

{~ l ) J 

V" 3 



7 = 1 J 7 = 1 i=0 V 7 1=0 7=1 W 7 



j°=l i=0 i=0 j=i 

J— 1 j J— i— 1 / , x „• J— 1 

= E|E ^v^E^'-D 

i=0 ' j=0 J ' i=l 

since E^^Tj 1 ^{i} — S/=i ^j; and we obtain X n CP(/i) as n — > oo. □ 

Second proof of Theorem [4] by Poisson approximation. By Lemmas [2] and HI we 
have 

n 

d(c(X n ), ^ CP(Bi(e fc ,e[ fc>n ]))J ^ d(Bi(e fc , Q[k, n ]), CP(Bi(e Jfc , g [Jfc)n ]))J. 

jt=i 

We prove that 

d(Bi(e,p),CP(Bi(e,p))) ^p 2 (E £ ) 2 (5.4) 

for all p G [0, 1] and for all non-negative integer-valued random variable e. By Barbour 
et al. [7J Corollary 10.L.1], we have 

d(Bi(e,p),CP(Bi(e,p))) < P(Bi(e,p) ^ l) 2 . 

Now 

oo oo 

P(Bi(e, p) > 1) = 1 - J> - p)< P(e = ^) = 5^ [1 - (1 - P) £ ] P(e = *) 

oo 
£=0 

hence we obtain (15. 4p . Applying (15.41) . we conclude 

n 

d{C(X n ), ^CPfB^efc^iM))) < 5^ m fc,i4,nl 

fc=i 

by (J3U). Clearly 

* CP (Bi(e A , £[*,„])) = CP I Bi ( £ fc> £[M) 
_1 \fe=i , 

hence, in order to prove the statement, it suffices to show 

n 

2_j BK £ fc, Q[k,n]) ~^ I 1 - 
fc=l 

We will check 

n 

P (Bi (e fc , Q[k,n]) = j) ~> ^0'} for a11 i e N - (5-6) 

fe=i 



(5.5) 



First note that by Taylor's formula, for all p G [0, 1] and all K, I G N, 



i=0 

where 



(i-p) x = E(f)(-i)V+^,/(p), 



i^,/(p)i<(f)p J - 



Hence 



p (Bi(e fc , g [k , n] ) =j) = 22[ . J e[ fcj „](i - e[*,n])'~ 3 p ( £ fc = £ ) 



J—j — l oo 



E E 6 \kri P ( £ k =Q + Rn,k,3 



i=0 ^=j+i • 
i=0 



where the sum is if j ^ J and 



~ /A 

-Rn,fc,j := > , f . 1 P[ k>n ] p { £ k = £)Rt-j,J-j{Q[k,n])- 

0—A \J ' 



Assumption (ii) implies ( 15. 3D again and we have 

E-j E — — m ^+^M = 7T E — T-E m w^[M^-y E —£T X 3+i = v{3} 

fc=l i=0 • y ' 8=0 ' fc=l J ' i=0 

as n — > oo for j — 1, . . . , J — 1. Moreover, (I5.3P implies 
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E i^i < E E y i,n] p (^ = o _ J 43 = E ^m^m -> ^ J > = °> 

hence we conclude (15.61). □ 



Next we study the case when the intensity measure /i of the limiting compound Poisson 
distribution CP(/x) may have unbounded support. 

Theorem 5 Lei (X n ) n( zz + be an inhomogeneous INAR(l) process. Assume that 

oo 

(i) Q n < 1 for all n6N, lim g n = 1, (1 - g n ) = oo, 

(ii) lim . \ = A,- G [0, oo) for all j G N suc/i £/iai i7ie limits 



J i=0 



exist. 



Then 



X n —> CP(/i) as n — ► oo. 



Proof. We follow the second proof of Theorem [4] by Poisson approximation. We have to 
show that fi is a finite measure on N and to check (15. 6p . First note that by Taylor's 
formula, for all p G [0, 1] and all K,I G N, 



E(f)(-i)V«(i-rt^E(f)(-i: 

o— n \ / i— n \ / 



i=0 v 7 i=0 

Hence for all / G N, 



P (Bi(e fc , Q [kM ) = j) = 22 ( • ; J i.njl 1 - ^IM^ P ( e fc = ^) 

£=j V/ 



£=7 XJ/ i=0 

E E jiii^-j-i)! 0IM p ( £fc = £ ) 

i=o £=j j y j / 



21 

(.-I)" 



7 z — ' z 

J i=0 



One can easily check that ( 15. 3ft holds for all j G N, and we obtain 

- 1 27 (-1 

limsup ^2 P ( Bi (^ 0[fc,n]) = j) ^ ~7 E — T 

n^oo , i ■ i 

fe=l 

In a similar way, for all J G N, 



7! * — ' Z! 

J i=0 



« 1 27-1 (-1)* 

liminf V" P (Bi(e fc , g [k ,n\) = j) > -r Y] 

n— >oo *— • 7! ' Z: 

k=l J i=0 

hence by the existence of the limits (15.71) we conclude (I5.6p . 
Finally, for all J G N, we have 

J J n n 



^2 Vki} = E i™o E P ( Bi ( £ fc, Q[k,n]) = j) = Hm ^2 P ( X ^ Bi ( £fc > ^ n l) ^ J ) 

j=l j=l k=l k=l 

n n 

^ lim > P (Bi( 

n. — >ca ' * x n — >no g * 



n— >oo f * n— >oo ■ 

fc=l fc=l 



using again (15. 5p . Consequently, YlJLi A*{j} ^ -^l < °°5 hence the measure /i is finite. □ 

Remark 3 A possible limit measure CP(/x) in Theorem [5] is a special compound Poisson 
measure, since its intensity measure \i has finite moments. Indeed, for all J, t G N, we 



have 

J n J 



Y.lU 1) •••(./ l)fi{j} = lim VV./U D---U l)P(Bi (e k , g [k , n] ) = j) 



n — >oo 

k=i j 



^ lim V V j(j - 1) • • • (j - £ + l)P(Bi (e k , g [k , n] ) = j) . 

n— >oo ^— ' ^— ' 

fe=l j=£ 

It is easy to check that for all p G [0, 1] and for all non-negative integer-valued random 
variable e we have 

oo 

j^jO'-l)---0'-' + l)P(Bi(e,p)=j)=p < Ee(e-l)...(e-< + l). 
Consequently, 

J n 

j(j - 1) • • • 0' - t+ !Mj} < lim Y] m k ,£Q[ k ,n] = Xe<oo 

j=l fc=l 

using again (15. 3p . 

Example 1 For n G N, let Qn = l-\ and P(e n = j) = -^j^, j G N, P(e n = 0) = l-± 
Then Ee n = oo for all n G N, thus inequality (15.41) is not enough to prove the compound 
Poisson convergence. Moreover, ^f"--'? = -4 for all j, n G N. The measure a on N 
defined by := -4, j G N, is finite and the infinite series Y^jLij^iJ} diverges. 

We prove that (X n ) neZ+ converges to CP(/i) in spite of the fact that assumption (ii) of 
Theorem O does not hold. We have, for n G N and z G C with \z\ < 1 and 2^0, 

TT . , 1 1 1^ z j l-zf, 4^ z* \ 1 (l-z)ln(l-z) 

= 1 + - > — r = 1 1 + > = 1 + 

n n ^ J \J + 1 ) n y j + 1 J nz 

which representation is valid on the whole D. By Lemma [T] we have 

T-r , (l-z)ln(*(l-z)) , 



fc=i 

Consider the functions F n : D — > D, n G N, defined by 



fnW :=n e " 



fc=l 

We have 

•fc 

■ J ""V ~'h l-f(l-z) 

as n — > oo. Since for the dilogarithm, see Abramowitz and Stegun [H Section 27.7] 



Li 2 (z) := \ — = - / — i i du, z <E D, 

~[ 3 Jo u 



holds, we have 



F n (z) -> exp <^ 



z 3 - 1 



as n — ■> oo 



for all z £ D. On the other hand, one can easily check that 
!*(£(!-*)) 



max 



as n — > oo for all z G -D with 1. (5-9) 



n(l- J(l-z)) 

Namely, for all z E D with 2^1, all n ^ 2 and all 1 ^ ^ n we have 





2 i 2afc 






l--(l-z 




H) 


-— i-N 


n 


n 







2a(n - 1' 



a: 



n- 



^ 1 - - < 1, 
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where a := 1 — Re z G (0, 2]. Moreover, 



— , for all ii G C with \u\ < 1. 



3=1 



Hence, for all z G D with z ^ 1, all n ^ 2, and all 1 ^ ^ n we conclude 



^ (|(1-^)) 



n. 



(i-5(i-*)) 



1 00 1 

n 7 



1-- l-z 

n 



-y-(i--) °'~ 1)/2 = ln £ - _> 

n ^-^ j \ nJ n x /l — r 



as n — > 00. An application of the inequality |e u — 1 — m|^|m| 2 valid for all w G C with 
|w| ^ 1/2 implies 



|F B (z)-F„(z)|< J] 



fe=i 



(l-*)ln(i(l-*)) 



n(l-±(l-*)) 



as moo 



for all z G D by (15. 8p and (15.91) . Thus we finished the proof. 

Open Problem. The above example shows that in Theorem [5] we do not exhaust 
the possible limiting compound Poisson distribution. We conjecture that every compound 
Poisson measure can appear as a limiting distribution of an inhomogeneous INAR(l) process. 

Theorem [3] can also be extended for the case of limiting compound Poisson distribution. 

Theorem 6 Let k n G N for all jigN, and {Cn,fc : 1 ^ k ^ k n , jigN} be non-negative 
integer-valued random variables with factorial moments 

m nAj := EC„, fc (Cn,fc - 1) • ■ ■ (Cn,k - j + 1), j e N. 
Moreover, let G [0, 1], 1 ^ k ^ k n , n G N. Assume that 



(i) S Pn fc m n,fc,j — *■ <\j £ [0) 00 ) / or a ^ J N suc/i t/iai t/ie /zmzts in (15.71) exist; 
fe=i 



(ii) Y. G°n,fcEC„, fc ) 2 -> 0; 
fc=i 

as n — > oo. Tnen 
m /aw as n — ► oo. 



* Bi(£(Cn,fc),Pn,fe) -> CP(/i) 



fc=l 



6 Appendix 



Lemma 3 If ak,bk E D , k — 1, . . . ,n, then 



fc=l fc=l 



< |ofc - fefcl- 



fc=i 



Proof. The statement follows from 

n n n /k—1 \ / n 

jj ctfe - jj &*, = 5^ ( n s ) ~ ^) ( n ^ 

fc=l fe=i fe=i \i=i / \j'=fc+i 
valid for arbitrary a^^eC, k = 1, . . . ,n. 

Lemma 4 If \i k , v k , k = 1, . . . , n, are probability measures on Z + then 

n 

d( * fi k , * u k ) < y^d(n k , v k ). 

\ k=l k=l / ' * 



k=l 



Proof. The inequality 



k=l k=l ' 



□ 



easily follows from the definition of the total variation distance, where f\ denotes product 
of measures. By Barbour et al. [7J Proposition A. 1.1], we have 



(n n \ n 

k=l k=l ' k=l 



and we obtain the statement. 



□ 



In the proofs we use extensively the following lemma about some summability methods 
defined by the sequence (g n )nEN of the offspring means. 

Lemma 5 Let (g n )nen be a sequence of real numbers such that g n 6 [0, 1) for all n e N, 

oo 

lim g n = 1, and (1 ~~ Qn) = 00 ■ Put 



n=l 



l n,j ■- 



(1 — Qj) jQ q\ for n,j,kEN with j ^ n. 
e=j+i 



Then a^- ^ a^- for all n,j, fcGN with j ^ n, 



(!) n 

max al '• — ► U as n — > oo, 



and /or an arbitrary sequence (x n )neN 0/ rea/ numbers with lim i„ = 1 6 I, 



Eai*o-^i ~ * — as n^oo for a// fc£N. 
k 

Proof. For each J e N, we have the inequality 

^ max aiH ^ maxaiH + max a£H ^ max(l — o,) + max (1 — o,-) [ [ 

hence letting n — > 00, we obtain 

^ lim sup max a^- ^ max(l — Oj), 

since 

n n n 

^ jQ o^ ^ exp < — (1 — Qe) > — > as n — >■ 00. 
t=j+i ^ i=j+i > 

Now letting J — > 00 we get 

^ lim sup max af^ ^ 0, 

and we conclude (16.11) . 

By the Toeplitz theorem, in order to prove (I6.2p . we have to show 



n— >oo 
n 



(k) 


= 




(*) 


1 


'n,j 


~ V 



lim > 1 

i=l 

n 

El ( k )\ ^ 
|a„j| < 00 



nGN - 



for all fceN. By the assumptions, 



n s n 

£=j+i ^ <=i+i 

and we obtain (16.41) . Next we prove (16.51) and (16. 6p for k — 1. We have 



Qe) r — * as n^oo 



n n n n / n 



z>s = e^ 1 - n « = e n ft-iift ) = 1 -lift -> 1 

i=i j=i £=j+i i=i \e=j+i e=j 



as n 



by (16.31) . and lim ^2 fl l = 1 a l so implies that sup \ a n j\ = SU P S a n i < 00 • 

j=% J ' nGN 3=1 nGN j=l 

we finished the proof of the statement of the lemma in case k = 1. 



The aim of the following discussion is to show (16. 5p and (16.61) for all k ^ 2. Observe 
that 



n / n 



i -u<s='E[ n «j-n«s =E (1 -* n a 

t=\ j=i \£=j+i t=j J j=i e=j+i 

n / k 



e(e $<!-*)' n 

j=l \ j=l \ / / £=j + l 

fc 



Ql 



i=l i=2 V 7 j=l 

n n 

where, by (I6.3p . ^ lim Yl Qt ^ li m Yl Ql = 0- Moreover, 



< 5^(1 - QjY^a^ <: ^(1 - ej) i_1 aJ2 ^ as n -> oo for all i > 2 

3=1 3=1 

by the lemma for k = 1 and by the assumption lim g n = 1. Consequently, we obtain 
(I63j) and hence (JEEJ for all Jfe > 2. □ 

Lemma 6 Lei e be a nonnegative integer-valued random variable with factorial moments 

m k := Ee(e - 1) • • • (e - k + 1), fceN, 

mo := 1, and with generating function H(z) = E(z £ ), defined for z E D. If < oo 
for some k G N i/ien 

fc-i 

#(z) = V - lV + R k (z) for all zeD, 

j=0 J ' 



where 

jfe! 



<Trk-l| fc forallzeD. 



Proof. By mj = £ ^ - 1) - (£ - j + 1) P(e = £), 

1=0 

k—l oo /• k— 1 \ 

ftw = aw -E>-iy=E»'-E (»- 1 ) J ) p ( e 



j=o J e=o v j=o 

and by Taylor's formula for the function z i— »■ z we get 

jfe! 



^-Ef-V 1 )- 

3=0 V/ 



^ l|z - l\ k sup |£(£ - 1) • • • (£ - k + 1)(1 + 9{z - l) f- k \ 



6>e[o,i] 



for all z e D. □ 
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